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We use the low density approximation within the kinetic theory approach to vacuum creation of an e~e + 
pair plasma (EPPP) in a strong laser field in order to investigate the dependence of the observed EPPP on the 
form and parameters of a single laser pulse. The EPPP distribution function is calculated for an arbitrary time 
dependence of the electric field in the multiphoton domain (adiabaticity parameter 7 3> 1). The dependence 
on the field strength, the form and spectrum of the field pulse is investigated on the basis both analytical and 
numerical methods. The obtained results can be useful for examining some observable secondary processes 
associated with the dynamical Schwinger effect. 
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1 Introduction 

The short wavelength domain of the laser radiation is very prospective for experimental observation and inves- 
tigation of the dynamical Schwinger effect |[]J[3). The operating X-ray laser facilities J4][5] are able to generate 
only very small electric fields in the focal spot with the amplitude Eq ~ 10~ 3 — 1CP 2 E c , where E c = m 2 /e 
is the Schwinger critical field strength. However, there are perspectives to achieve the Schwinger limit in this 
domain in the nearest future (e.g., ]6]|2l). 

In the present work we investigate the residual electron-positron pair plasma (EPPP) which is generated from 
vacuum after the cessation of a high-power laser pulse in the focal spot. In contrast to the well-known works 
(HE), our approach is based on the analysis of a kinetic equation (KE) 110111 11 that is an exact nonperturbative 
consequence of the basic equations of motion of QED for the case of a linearly polarized time dependent electric 
field ("laser field"). We limit ourselves to the domain of subcritical fields Eq -C E c and to the case of rather 
short wavelength A when the adiabacity parameter 7 = 2ttE c X c / EqX is large, 7 ^> 1. This corresponds to the 
multiphoton mechanism of the EPPP excitation (A c = 1/mis the Compton wavelength). 

In problems of vacuum particle production, one usually studies the characteristics of the out-state, when the 
external field is switched off. This is due to the needs of an experiment aimed at finding free real particles. 
Recently, however, there is a growing interest in the properties of the intermediate states (mid-states), character- 
ized by quasiparticle interactions in the presence of an external field. This is motivated by the likely role of the 
elementary reactions of the quasiparticles in the field and their possible experimental manifestations |[T1 [T2HT4l . 

The most complete description of the pair production process for both the mid- and the out-states is given by 
the momentum distribution function. Depending on the shape of the external field the vacuum excitations can 
form very complex and diverse momentum distributions which are essentially nonequilibrium ones. The com- 
plex form of the pair momentum spectrum should significantly affect the cross sections of elementary reactions 
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and the development of avalanche-like electromagnetic cascades 11311141 . It should also have some observable 
manifestations. 

The most striking difference between the nonequilibrium spectrum of produced particles and the equilibrium 
one is a significant non-monotonicity of /(p). This deformation is much more pronounced in the mid-state 
which is important for elementary reactions and cascades. This feature has been noticed in the works JTJQ3] 
for the case of a periodic field, and was later studied more in detail for a Gaussian laser pulse in the work lfl6l . 
The representative oscillations of the momentum distribution with a scale set by the laser frequency received an 
explanation on the mathematical basis for the Stokes phenomenon in the theory of differential equations ifTTI . 
Other manifestations of such effects were studied in II18II19I . 

It is important to note that these remarkable features of the momentum distribution function are not captured 
by WKB-like approaches (8l|9). On the other hand, the only known exact solution for a homogeneous field in 
the form of a Sauter bell generates a smooth spectrum j|20l . Apparently, for the occurrence of nonmonotonic 
deformations of /(p) different scales must be present in the spectrum of the external field. The actual laser field 
is required in the spectrum of two different time scales, so that the question of the form /(p) is essential for the 
experiment. 

In the framework of using limitations, the following main results were obtained by both numerical and analyt- 
ical approaches: 

• the distribution function of the residual EPPP has a complicated quasiperiodical character in both the longitu- 
dinal (pii) and the transversal (p±) momentum components resulting in a cellular structure in the momentum 
space; 

• the profile function for switching-on (and switching-off) the laser signal has a considerable influence on the 
magnitude and structure of the momentum distribution of the residual EPPP; 

• together with the usual multiphoton mechanism of the EPPP excitations the photon cluster mechanism of 
absorption of the laser field energy plays an essential role due to the relativistic nature of this system. 

In Sect. [2] we explain the problem and give an analytic solution of the KE in the domain 7> 1 and Eq <C E c . 
Such analytical estimations are important due to the complications encountered in the numerical calculations. 
In Sect. [3] we present some numerical solutions of the KE in this domain. The basic results are summarized in 
Sect, g] 



2 Low density limit in the multiphoton domain (7> 1,£J < E c ) 

For the description of the EPPP vacuum creation we use the exact nonperturbative KE obtained in the work [10] 
for an arbitrary time dependent and spatially homogeneous electric field with linear polarization for which we 
use the Hamiltonian gauge A»(t) = (0, 0, 0, A(t)) and the field strength E(t) = -A(t), 

t 

/(p,i) = ^A(p,i) f dt'X(p,t')[l-2f(p,t')} coB0(t,O, (D 

to 

where 

A(p,i) = eE(t)e^/uj 2 (p,t) (2) 



is the amplitude of the EPPP excitations, w(p, t) = y £j_(p) + (p\\ — eA(t)) 2 with the transverse energy e± = 
(m 2 + p 2 ^) 1 ! 2 and the high frequency phase is 



0(t,t') = 2 J druj(p,T) . (3) 

The distribution function in the quasiparticle representation f(p,t) =< in | a + (p, t)a(p, t) in > is defined 
with the in-vacuum state | in >. For the generalization of the KE (fl} to arbitrary electric field polarization 

see nsjEUEa. 
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In the low-density limit /< 1, the KE (Q]) leads to the formal solution 



/iow(p,i) 



dt'X(p,f)e i9 ^ 



(4) 



Here, the external field is switched on in the infinite past, to — > — oo. The approximation implies that the 
external field is rather weak, E <C E c . Below we will restrict ourselves to the analysis of Eq. (|4j in the limit 
t — > oo, which defines the momentum distribution of the residual EPPP, 



/out(p) = lim /iow(p,i) = 7 



dt\(p, t)e 



iO(t) 



(5) 



Here the representation of the phase ((3} via antiderivatives has been used, 9(t, f) = 0(f) — 9(f). 

Hereafter we will consider the domain of the multiphoton mechanism of EPPP vacuum creation, where the 
adiabaticity parameter is large ||8], 7 > 1. A perturbation theory w.r.t. the small parameter 1/7 < 1 can be 
constructed here. 

As a first step let us consider the function / out (p) © in the leading approximation, when 



w(p,t) 
X(p,t) 



w 



(p) = Vp 2 



X (p)E(t), X (p) = ee±/w$(p) 



(6) 



and 9(f) = 2u>o ■ t. This means that the contribution of the high frequency harmonics of the external field is 
omitted in Eq. ©. It leads to the simple result 



/out(p) 



1 



A; 



E(u = 2u ) 



(7) 



where E(tu) is the Fourier transform of the field strength E(f). In this approximation the vacuum EPPP produc- 
tion takes place when for the frequency of the one photon e~e + pair creation process wi 7 holds wi 7 = 2wo- This 
mechanism is exclusive here. Its intensity is regularized by the presence of the frequency wi 7 in the spectrum of 
an external field. 

In order to facilitate EPPP creation, one can take into account the multiphoton mechanisms of EPPP produc- 
tion. These processes are described by the high frequency multiplier in Eq. (0. Let us use the non-perturbative 
method of photon counting assuming that the electric field is periodical, A(f) = A(t + 2%/v), where v = 2ir/T 
is the angular frequency and A(t) = A(—t), in order to provide the property of being an even function of the 
quasiparticle energy, uj(p, t) = uj(p, —f). The corresponding Fourier transform 



co(p, t) = 0„ cos nut 



(8) 



leads to the decomposition of the phase in Eq. (f5J 

9(f) = 2fl t + 

a n smnvt, 



(9) 



where a n = 2Vt n /vn. In the case of a harmonic external field 80 is the renormalized frequency J9]. Let us 
employ now the leading approximation A(p, t) —> Xo(p) E(f) (O and a non-perturbative decomposition based 
on the well known formula 



exp (m sin </>) = J n (a)e 



irufi 



(10) 



n— — 00 



where J n (a) is the Bessel function. 
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Let us now consider the integral in Eq. (0 and perform the substitutions d8ll-([T0li, leading to 

oo 

J = J dtE{t)e im 

— DO 
OO 

= j dtE(t)e 2mot (j (2n ) + I] E ^KK (fe "" W + (k n -> -k n )\ . (11) 

The first term with Jq(2£Iq) describes the direct vacuum excitations at the frequency uj = 2Hq and corresponds 
to Eq. (0. If we are interested in excitations at lower frequencies, this contribution can be omitted. 

A new feature of the multiphoton processes here is the cluster character of the energy absorption from the 
photon reservoir of the external field. The photon cluster of n th order consists of n identical photons with the 
carrier frequency v and the total energy nv. As it can be seen from Eq. (fTTl i. it is possible that k of the ?i-photon 
clusters are absorbed simultaneously. The ordinary multiphoton process corresponds to the simplest "cluster" of 
the order k = 1. The appearance of photon clusters in the multiphoton processes in Eq. (fTTT i is a consequence 
of the nonlinear field dependence of the quasiparticle energy w(p, t). The probability for the generation of an 
n-photon cluster is defined by the amplitude a n in the decomposition (O (the argument of the Bessel function 
Jk{o-n) in Eq. ( fTTT i) while the probability of the simultaneous fc-cluster absorption corresponds to the Bessel 
function of the order k. 

We rewrite the integral ([Til in the following approximation: 

oo 

/ dtE(t) c 2lUot J [ I Jk„{a n ) cxp I i ^ (nk n )vi J + (k n -> -k n ) \ , (12) 

"=1 k n = l ^ ^ n=l ' ' 

where n max is the maximal photon number in the cluster, J = </(n max — )• oo). Let us perform here the resum- 
mation procedure. In Eq. (fT2l we conserve the index k n for notation of the k n identical photon clusters where 
each cluster contains n photons. Then the total photon number in this photon set will be equal to Nk = nk n . Let 
us replace now the summation over n by a summation over Nk. The integral (fT2l can then be rewritten as 

JV m „ 

J = 2ttJ2 e (u = 2^o - v J2 l N k)) II J k( a N k /k) ■ (13) 

fc=l N k = l JV fc = l 

Here denotes an arbitrary integer multiple of k and N max ~ n max in Eq. ( fT2l . 
Substituting Eq. fl3i into Eq. ((5]), we obtain the generalization of Eq. (0 



/out(p) =7T 



2V, 

2i2' 



fc=l iV fe =l iV fc = l 



(14) 



The summation is developed here by the order of the photon cluster. Below we will restrict the sum to the usual 
multiphoton processes (k = 1). 

Thus, the multiphoton mechanism leads to a lowering of the EPPP excitation frequency. The maximal effect 
is reached on the carrier frequency v of the external field, i.e. for the condition 



2fi - v J2 [Nk] = v- (15) 



N k =l 



Let v <C m. Then the hypothesis of the phase randomization can be used, 

~< E{Nv)E*{N'v) \E{Nis)\ 2 8 NN , , (16) 
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where N, N' are integers and -< • • • ^Au denotes the averaging procedure over a small frequency interval Aui. 
The approximation (fT6l l applied to Eq. ( TBI leads to the result 



k=l 



N k = l 



In the case of the monochromatic external field Eq. (TPTT i exhibits the accumulation effect. Using the formula 
6 2 (lu) = (T/2n) ■ 6(u) fovT > 2%/u, one can obtain from Eq. ([17) the EPPP production rate 

Wp) = Ag E II •>* - f £ ™ - ") ' ( 18 > 

fe=lAT fe =l JV fc = l 

Then the condition (TT3T > defines the photon number required for the breakdown of the energy gap with the 
renormalized mass m rcn = f2rj(p = 0). Since the renormalized frequency f2o(p) is a monotonously increasing 
function of p = |p| with min fio(p) = wiren> this photon number will be minimal for the particles created at rest, 
Nq = N max (p = 0). In order to investigate the momentum distribution of the EPPP it is necessary to consider a 
larger photon number, N max > Nq and iV max — > oo forp — > oo. 

The simplest situation corresponds to the usual multiphoton processes when the photon cluster mechanism is 
a minor, k = 1. The Eq. ( TT~8T > in this case has the form 



/out(p)= 5 [2fi (p) - K^max + 1)] JJ J lM ' ^ 

n=N 

The set of the Bessel function arguments a n corresponds here to the spectral composition of the multiphoton 
process in Eqs. © and (|9). 

Already on the qualitative level it can be seen from Eqs. (fT8T l and ( fT9] l that the momentum distribution of the 
residual EPPP has a cellular structure in momentum space stipulated by the presence of the Bessel function. It 
leads to some quasi-periodic behavior of / ou t(p) w.r.t. the components p» andpj_. 

Eqs. (TT~8T > and $19[ allow to perform various estimations of EPPP characteristics (particle number density, 
momentum distributions) for different parameters of the laser radiation in the domain 7 ^> 1. 

As an example we find the estimate for the total intensity of EPPP production 

/ = 4 /(2^ /out(p) (20) 

for v — m and Eq = 0.1 E c . It corresponds to the two-photon mechanism of EPPP creation. From Eq. ( fT9] l it 
follows then that the rate of EPPP production in the volume A 3 per period of the field action (Eq /E c = £) is 

I = 1-S 2 J? (C 2 /4)m 4 w 10~ 9 /(A 3 r) . (21) 

Z7T 



3 Numerical calculations 

We explore field pulses with two types of shape: the radio pulse (N periods of the harmonic field) 

E(t) = E sin ut, < t < N 2tt/u, N - integer, (22) 
and the Gaussian pulse 

E(t) = Eq exp (-t 2 /2r 2 ) cos vt, vt = a, -10r < t < lOr . (23) 

Figures [T]|3]|4] show the complicated behaviour of the momentum distribution / ou t(p) f° r tne periodical field 
(1221 in dependence on the field parameters. This distribution turns out to be asymmetric w.r.t. the interchange 

Copyright line will be provided by the publisher 



6 



D. B. Blaschke et al: Influence of Laser Pulse Parameters on the Properties of e e + Plasmas 





Fig. 3 Increasing complexity of /out(p) in the field ( 122b with increasing TV at Eo — 0.01 E c and 7 = 24. Left panel: N=l; 
right panel: N=2. 



P| I ~P\\ (Fig- 03- Increasing the number TV of periods in the pulse (|22] > leads to a complication of the 
momentum distribution (Fig. [3) and to the generation of a cellular structure for large N ^> 1 (Fig. |4j. These 
details are lost when using the WKB methods [H|9l ■ 

Analogous alterations can be observed in the case of the Gaussian pulse d23l l. For example, Fig. |2]shows the 
distribution function for a = 5 and 7 = 10. Fig.[5]illustrates the strong dependence of n out on the field strength. 
A comparison of the EPPP number density for the pulses (l22l i and d23l l shows that the periodical field can produce 
a significantly higher output of observable pairs. The obtained estimates for n out are notably different also from 
the well-known results of works based on the WKB approach. 
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Fig. 5 Dependence of n out on the field strength for a Fig. 6 Comparison of n ou t for different pulse shapes and 

Gaussian pulse i23\ for two values of 7. different theoretical approaches. 
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4 Conclusions 

In the present work we have done first steps in the investigation of the residual EPPP based on a kinetic theory 
foundation for strong nonperturbative QED. We restricted ourselves here to the multiphoton domain 7 3> 1 
and the low-density approximation. We have shown that the momentum distribution of the residual EPPP is a 
strong nonequilibrium one and has a complicated cellular structure in the momentum space depending on the 
characteristics of the laser radiation. As a result of the numerical analysis it was observed that the residual 
EPPP density depends strongly on the pulse shape for switching on (and switching off) the laser field. Finally, 
a theoretical approach was developed for estimating different characteristics of the residual EPPP. This method 
allows to take into account both the usual multiphoton processes of EPPP excitation and also the new photon 
cluster mechanism. 
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